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The  Energy  Release  Rate  (ERR)  for  the  quasi-static  problem  of  a 
semi-infinite  mode  I  crack  propagating  through  an  inhomogeneous  isotropic 
linearly  viscoelastic  body  is  examined.  The  shear  modulus  is  assumed  to 
have  a  power-law  dependence  on  depth  from  the  plane  of  the  crack  and  a  very 
general  behavior  in  time.  A  Barenblatt  type  failure  zone  is  introduced  in 
order  to  cancel  the  singular  stress  and  a  formula  for  the  ERR  is  derived 
which  explicitly  displays  the  combined  influences  of  material 
viscoelasticity  and  inhomogeneity.  The  ERR  is  calculated  for  both 
power-law  material  and  the  standard  linear  solid  and  the  qualitative 
features  of  the  ERR  are  presented  along  with  numerical  illustrations. 


1.  Introduction. 


In  an  earlier  paper,  [1],  the  quasi-static  problem  of  a  semi-infinite 
mode  I  crack  propagating  through  an  infinite  inhomogeneous  isotropic 
linearly  viscoelastic  body  was  investigated.  The  primary  focus  of  that 
study  was  the  description  of  the  stress  and  displacement  fields  in  a 
viscoelastic  material  characterized  by  a  shear  modulus  exhibiting  a  very 
general  behavior  in  time  and  a  power-law  spatial  dependence.  The  principal 
objective  of  this  paper  is  to  calculate  the  Energy  Release  Rate  (ERR)  for 
the  crack  problem  studied  in  [1]  and  to  illustrate  its  dependence  on  the 
combined  influences  of  material  viscoelasticity  and  inhomogeneity. 

In  this  work,  the  viscoelastic  body  in  which  the  crack  propagates  is 
assumed  to  have  a  shear  modulus  of  the  form 


k*  (  C ,  y )  =  Wcg(t)(y/yc> 


(1) 


where  y  denotes  distance  as  measured  from  the  plane  of  the  crack,  y 
is  a  characteristic  depth,  0  <  Y  <  1,  and  initially  the  only  restrictions 
placed  on  g  are  that  it  be  a  positive,  continuous,  decreasing,  and  convex 
function  of  time.  Subsequently,  for  illustrative  purposes  the  specific 
forms  of  g  to  be  considered  will  correspond  to  the  standard  linear  solid 
and  a  power-law  material. 

The  assumption  of  a  simple  power-law  behavior  in  space  for  the  shear 
modulus  introduces  certain  unphysical  features  into  the  model.  However,  as 
pointed  out  in  [1]  and  [2]  and  references  cited  therein,  the  consequences 
of  the  pure  power-law  model  do  not  preclude  its  applicability  to  certain 
physical  scenarios.  Apart  from  the  suitability  of  the  shear  modulus 
considered  in  this  paper  to  any  specific  model,  an  important  aspect  of  this 
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work  is  the  demonstration  that  for  such  a  modulus,  a  tractable  boundary 
value  problem  results,  from  which  can  be  obtained  meaningful  average 
quantities  associated  with  the  energetics  of  the  crack,  specifically,  the 
ERR.  Furthermore,  it  shall  be  shown  by  a  simple  perturbation  argument  that 
for  materials  which  experience  a  significant  softening  near  the  plane  of 
the  crack,  a  notion  to  be  made  precise  later,  the  simpler  assumption  that 
the  shear  modulus  vanishes  in  the  plane  of  the  crack  may  still  provide  a 
model  for  making  acceptable  engineering  predictions. 

The  Stress  Intensity  Factor  (SIF)  and  the  ERR  are  familiar  and 
important  notions  in  the  development  of  a  fracture  criterion.  For  a  large 
class  of  crack  problems  in  homogeneous,  linear  elastic  materials,  a  simple 
relationship  exists  between  the  two  quantities.  However,  in  the  case  of 
crack  propagation  in  an  inhomogeneous  viscoelastic  material  as  considered 
here,  the  usual  notion  of  the  SIF  is  no  longer  valid  since  the  order  of  the 
dominant  singularity  in  the  stress  field  near  the  crack  tip  is  a  function 
of  the  spatial  inhomogeneity  exponent.  However,  the  notion  of  the  ERR 
still  provides  a  meaningful  fracture  criterion.  Furthermore,  by  adopting 
the  techniques  and  loading  assumptions  utilized  in  [3],  the  ERR  is  easily 
calculated  and  illustrates  vividly  the  spatial  and  viscous  effects. 

The  boundary  value  problem  considered  here  is  that  of  a  mode  I 
semi-infinite  crack  propagating  with  constant  speed  in  a  nonhomogeneous 
isotropic  viscoelastic  media  in  a  state  of  plane  strain.  If  the  crack  is 
assumed  to  propagate  in  a  plane  about  which  the  spatial  properties  of  the 
body  are  symmetric  and  along  the  Xj-axis  with  speed  v,  driven  by  loads 
f(Xj-vt)  which  follow  it,  then  the  specific  boundary  value  problem  to  be 


solved  is 
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Here  o  e  ,  and  u.^  are  the  viscoelastic  stress,  strain  and 
displacement  fields,  v  is  Poisson's  ratio  (assumed  to  be  constant),  p 

is  the  shear  modulus,  <5^  is  the  Kronecker  delta,  and  p*de  denotes  the 
Riemann-Stieltzes  convolution. 


P*de  =  IL  W(t-T)  de(T). 


After  adopting  the  Galilean  variable  x^  =  x-vt  and  the  change  of 
variables  y^  =  y,  =  o^,  o^  =  a  =  u^,  etc.,  it  is  useful  to 

record  for  future  reference  the  following  results  derived  in  [1]: 


a 

yy 


(P,0) 


"Cq t ipvg(-vp) ] 


i-y 


Uy(P  ,0) 


p/CY+2)  cos(Y"/2) 
q(l-v)  I  sin(qir/2)  y^ 


q  =  [ (  1+Y ) ( 1-YV /( l-v ) ) ] 1/2 

I  =  2Y(Y+2)  B(I±a±l,  Iza±3) 


(2) 
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COS(Ytt/2)  r° 
,  x(l-Y)/2 


o  (s ) 
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j  s  j ( 1 - Y ) / 2  ds 


x>0 . 


(3) 


From  (3)  it  is  easy  to  show  that  near  the  crack  tip,  as  x+0,  the 
singular  terra  in  the  asymptotic  expansion  of  the  normal  stress  is  given  by 


(x) 


K(Y) 

x(l-Y)/2 


where 


K(Y  )  = 

In  the  above  equations  the 
defined  by 


-cos(tty / 2)  r ^  ° yy^s ^ _ 

7T  J-oo  (l+Y)/2 


ds , 


Fourier  transform  f  of  a  function 


f 


(4) 

is 


f(p) 


r00 

=  /_«  f(x) 


eixp 


dx , 


with  inverse 

fv(p)  =  /_„  f(x)  e-lxp  dx  , 

and  f+(x)(f  (x))  denote  the  restriction  of  f(x)  to  x>0(x<0). 

In  the  next  section  of  this  paper  a  simple  expression  for  the  ERR  will 
be  derived  for  a  general  form  of  the  time  dependence  in  the  shear  modulus. 
The  derivation  of  this  formula  will  parallel,  within  a  more  general 
setting,  the  calculation  of  the  ERR  carried  out  in  [3].  In  the  final 
section  of  the  paper  the  ERR  is  calculated  for  some  specific  models  and  the 
qualitative  features  of  the  ERR  are  discussed.  The  paper  concludes  by 
illustrating,  by  means  of  a  perturbation  argument,  the  effect  upon  the  ERR 
of  modifying  the  spatial  inhomogeneity  in  (1). 


Section  2.  The  Calculation  of  the  Energy  Release  Rate. 


For  fracture  in  linearly  elastic  material,  the  ERR,  hereafter  denoted 
by  G,  is  usually  defined  as  the  difference  between  the  total  power  input 
to  the  body  by  external  forces,  P,  and  the  rate  of  change  of  total 

•  c 

internal  strain  energy,  E,  and  kinetic  energy,  K.  G  can  then  be 
calculated  from  knowledge  of  the  singular  asymptotic  stress  field  at  the 

crack  tip.  For  a  viscoelastic  material  the  difference  P-E-K  now  involves 
a  term  representing  the  energy  dissipated  through  viscous  effects  and  the 
quantity  G  which  depends  upon  the  entire  history  of  the  singular 
asymptotic  stress  and  strain  fields  at  the  crack  tip.  Because  it  is 
difficult  to  calculate  G  from  the  singular  fields,  a  Barenblatt  type 
failure  zone  is  introduced  behind  the  crack  tip  in  order  to  cancel  the 
singular  asymptotic  fields  in  front  of  the  crack.  More  precisely,  the 
boundary  value  problem  presented  above  is  modified  by  assuming  that  two 
loads  are  acting  on  the  crack  faces:  the  applied  external  tractions 

o  (x,0)  =  f(x),  and  now  denoted  o  (x),  and  cohesive  failure  stresses 
y  y  s 

af(x)  acting  in  a  failure  zone  of  length  af  immediately  behind  the 

crack  tip.  The  only  assumptions  about  a^(x)  are  that  af  is  small 

relative  to  some  length  scale  a  associated  with  a~(x )  and  that  if 

e  e 

K  (Y)  and  K- (Y )  are  the  coefficients  given  by  (4)  corresponding  to 

a  and  o  respectively,  then  K  (Y)  +  KC(Y)  =  0.  Hence  the  effect 
e  t  e  f 

of  the  failure  zone  is  to  cancel  the  singular  stresses  ahead  of  the  crack 
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tip  and  thereby  produce  a  cusp  shaped  crack  profile  behind  the  crack  tip, 
What  results  is  that  for  the  steady-state  problem  considered  here  G  is 
given  by 


G  =  -/°  af(x)  u  (x,0)  dx 
a  L  y  y A 


(5) 


where  now  u  (x,0)  is  the  crack  face  displacement  corresponding  to  the 


combined  loading  o  +  a 

e  f 


A  significant  simplification  in  the  calculation  of  (5)  occurs  for 
special  cases  of  0  and  o  .  Specif icially ,  the  external  load  o~ 


and  failure  zone  stress,  o  will  be  assumed  to  have  the  fo 


rms 


x/a 


x/ a. 


ae^x^  =  Lee  *  af(x)  =  “Lfe  ,  <  x  <  0.  (6) 


For  af/ag  «  1,  the  essential  features  of  the  Barenblatt  model  are 
satisfied  by  the  assumptions  (6),  namely,  a  set  of  cohesive  stresses  and 
associated  length  scale  af  and  a  length  scale  ag  associated  with  the 

applied  load  such  that  cancels  the  singular  stresses  produced  by 

°e*  °f  an<i  °e  are  given  by  (6),  (5)  is  replaced  by 


G  -  „  o  (x)  u  (x,0)  dx. 

L  y  > x 

A  general  method  for  calculating  (7)  begins  with  (2).  If 


T(p)  =  “Cq [ipvg(-vp) ]  p 


1-Y 


,  then  (2)  may  be  written  as 


(7) 


a  (p,0)  =  T(p)  u^( p,0)  -  a  (p,0). 


(8) 
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With  the  introduction  of  the  Barenblatt  failure  zone  u^  represents  the 
displacements  corresponding  to  the  combined  external  and  cohesive  loads, 

a  =  a  =  a  +  a  and  o*  =  o+  is  the  resultant  stress  ahead  of  the 
yy  e  f  yy 

+ 

cracks.  The  solution  of  (8)  is  given  in  [1].  Specifically,  if  F  (z), 
functions  analytic  in  the  upper  and  lower  complex  half-plane,  are  defined 
by 


F  (z) 


X~(z) 

2*i 


-q  (Q 

x+(x) 


dT 

T-z 


(9) 


with 


X+(p)  =  T(p)  X  (p). 


(10) 


then 


and 


0+(p)  =  lim  F+(z) 
lm(z)*0 


u  (p)  =  lim  _  F  (z). 

^  lra(z)*0 

Furthermore,  X+(z)  has  a  particularly  simple  analytic  extension  to  the 
upper  half-plane  given  by 


X+(x) 


=  c  / g(0) 


( 1  -T  )  /  2 


(11) 


c  a  constant. 

Applying  the  Parseval  formula  for  the  Fourier  transform,  (7)  may  be 


written  as 
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G  "  oo  5f(p)  uy>x(p)  dp- 


From  (6)  it  follows  that 


~afLf 


-af  Lf 

,  ®f(p)  =  u-  . .  . - 7 

f  2it(l-ia,p) 


03) 


with  similar  expressions  for  a  and  o  .  Since  u  has  an  analytic 

e  e  y  J 

extension  to  the  lower  half-plane,  namely  F~(z),  and  -ipu  ■  u 

y  y,x* 

it  follows  that  u  can  also  be  extended  analytically  into  the  lower 

y 

half-plane  and  from  (13)  it  is  clear  that  o has  a  meromorphic  extension 
to  the  lower  half— plane  with  a  simple  pole  at  — i/a^.  Hence  the  integral 
in  (12)  may  be  evaluated  using  residues  from  which  it  follows  that 


G  =  -  —  F  (-i/a,), 
af  t 


(14) 


To  evaluate  F  (-i/a^),  begin  by  writing  F+(z)  =  F+(z)  +  F*(z) 
where  from  (9) 


F+(z)  =  r  qe(T)  dr_ 

6  2ui  ^  X+(t) 


(15) 


The  analogous  formulas  to  (11)  for  o  and  a  show  that  a~  has  a 

e  e  e 


simple  pole  at  i/ag  and  since  X  (z)  is  analytic  in  the  upper  half¬ 
plane  (15)  may  be  evaluated  using  residues  to  give 
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a  L 


F+(z)  =  -X+(z)  {— - i _ J 

P  V/Vkl  + i  a  7  V  4-  4.  ' 


1+ia  z'  x  +  ,  , 

e  X  (z)  X  (i/a  ) 

e 


and  so 


Fg(p)  =  P)  U  -  + 


X+(p) 


x  (i/a  ) 
e 


Similarly  it  follows  that 


Now  one  may  write 


Ff(P>  =  -MP)  (l  -  — X--(P)  j 


X  (i/af ) 


°.(P)  =  a*(p)  +  o"(p) 
^  “  6 


=  F^(p)  +  ae(p) 


=  ag( p) 


X+(p) 

X"(i/a  ) 
e 


and  similarly 


°f(p)  =  of(p)  + 


X+(P) 


X  (i/af ) 


It  now  follows  that 


°(p)  =  °e(p)  +  °f(p) 


X+(p) 


a  L 
e  e 


afLf 


X+(i/ae)(l+iaep)  X+(i/af )(l+iafp) 


06) 
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From  (4)  and  the  requirement  that  Kg(Y)  +  Kf(Y)  =  0  there 


( 1-Y )/2  _  ( 1-Y )/2  . 

3e  Le  "  af  Lf» 


or  equivalently, 


X+(i/ae)  X+(i/af) 


Thus  (16)  may  be  written  as 


S(p)  -  (araf)  L*  ^  . 

( l+iaop)(l+iafp)  X+(i/a  ) 


Since  F  (p)  =  uy(p)  =  o(p)/T(p),  combining  (14),  (17),  and 
follows  that 


G  = 


(af"ae)  Le  X_(-i/af>  X+(i/af) 


(a  +a  )  2 
r  e 


[x+(i/ae)j 


Using  (10)  and  (11)  there  finally  results 


„  -  (Vaf)  Le  (1-Y ) 

<Vaf>  2C0  a 


results 


(17) 


(18) 


(18)  it 


(19) 


(20) 


vg(iv/af ) 


1 1 


It  should  be  noted  that  in  general,  beginning  with  a  relationship  as  in 

(8),  G  will  be  given  by  (19)  where  X  (z)  are  solutions  to  the  equation 
(10).  As  shown  in  [3],  the  above  derivation  can  be  modified  to  produce  an 
analog  of  (20)  for  more  general  loads  of  the  form 


®e<x) 


°f(x) 


tx/ a 

e 

e 


-  -h  Jo 


tx/ a 
e 


dhe(t) 

f  dhf(t) 


where  he(t)  and  h^(t)  are  arbitrary  signed  (not  necessarily  positive) 
measures  restricted  only  to  the  extent  that  the  required  integrals 
converge.  In  particular,  the  above  derivation  for  G  corresponds  to 
dhe(t)  =  dhj(t)  =  6(t-l),  the  Dirac  measure  concentrated  at  t=l.  However, 
for  the  more  general  loads,  the  analog  of  (20)  presents  more  numerical 
complications  in  the  calculation  of  G. 

Section  3.  Examples  and  Results. 

Two  special  cases  of  g(t)  corresponding  to  the  standard  linear  solid, 

g(t)  =  1  +  n  exp(t/x),  and  a  power-law  material,  g(t)  =  1  +  (t/x)”n, 
where  t  is  a  characteristic  relaxation  time  will  now  be  considered.  If 

the  Carson  transform  g(t)  of  the  function  g(t)  is  defined  by 


g(s)  =  g(0)  +  / "  e  ts  dg(t). 


and  the  following  nondimensional  parameters  are  introduced 


=  af/ae>  6  =  a  =  Tv/af, 


c  e’ 
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then  after  recalling  the  definition  of  and  noting  that 

g(ivp)  *  ivpg(-vp)  (20)  becomes 


G  = 


(l-e)  aeLe  (2<S)T  l 


(1+0  2p  E(Y)  . 

c  g(v/a  ) 


(21) 


with 


E(Y)  = 


[T (y-t-2)  ]  cos(Ytt(2) 


q(l-v)  sin(qu/2)  r(I±S±I) 


For  the  standard  linear  solid 


g(v/a  ) 


l+a( 1+n ) 

1-KX 


and  for  power-law  material 


(22) 


g(v/af)  =  l+o  r(l-n), 


(23) 


The  numerical  calculation  of  the  ERR  from  (21)  is  a  simple  matter,  and 
certain  qualitative  features  of  G  are  immediate.  The  inhomogeneous 

v 

elastic  effects  in  (21)  are  contained  in  the  term  (2<5)  /E(Y)  while 

l/g(v/af)  represents  the  viscous  effects.  The  purely  elastic  problem  is 
recovered  when  the  crack  speed  v=0  and  the  homogeneous  problem 
corresponds  to  Y-0.  Consequently,  the  qualitative  and  quantitative 
effects  of  adding  material  viscoelasticity  and  inhomogeneity  to  the  basic 
elastic  model  are  easily  identified  and  studied. 

From  (21)  it  is  easy  to  determine  an  "effective  depth",  y  ^ ,  such 
that  the  ERR  predicted  by  the  inhomogeneous  material  model  equals  the  ERR 


13 


predicted  by  a  homogeneous  model  with  shear  modulus  equal  to  the 
inhomogeneous  modulus  evaluated  at  y  =  y^.  For  a  homogeneous  material 
with  modulus  given  by  the  inhomogeneous  model  evaluated  at  an  effective 
depth  y^  ,  i.e.,  for  which 


u(t)  =  g(0  =  g(t), 


the  ERR  is  given  by  (21)  with  Y=0  and  P  substituted  for  p  .  Hence 

e  c 


G  = 


1-e 
1+e  • 


a  L2 
e  e 


2ue(l-v)  g(v/af) 


(24) 


Equating  G  in  (21)  and  (24)  there  results 


yl  2  1  l-u  1 


v 

To  illustrate  the  result,  (y^/a^)  is  graphed  in  Figure  1  when  Poisson’s 
ratio  v=,3. 

It  is  apparent  that  G  is  a  monotonically  increasing  f  uncion  of  <5  • 
This  behavior  is  explained  by  noting  that  as  6  increases,  so  does  the 
distance  from  the  plane  of  the  crack  at  which  the  material  is  as  rigid  as 
the  homogeneous  material  with  modulus  d(t)  =  g(t).  Thus  with 

increasing  6  there  is  a  larger  region  of  ’soft’  material  about  the  plane 
of  the  crack.  In  previous  studies  of  crack  problems  in  nonhomogeneous 
elastic  materials  (c.f,  [4]  and  the  accompanying  references)  it  was 
observed  that  with  a  decrease  in  rigidity,  there  is  an  attendant  increase 
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in  crack  instability.  Since  G  has  the  interpretation  of  power  available 
to  the  crack  tip  for  propagating  the  crack,  a  monotone  increasing 
dependence  on  6  is  in  qualitative  agreement  with  these  previous  studies. 

It  follows  immediately  from  (22)  and  (23)  that  the  ERR  is  a  monotone 
decreasing  function  of  a  and  hence  of  crack  speed.  Regardless  of  which 
model  is  used,  the  qualitative  differences  in  G  due  to  viscous  effects  are 
similar.  For  this  reason  the  normalized  ERR,  g,  defined  by 


a  L2 
e  e 

2u 


g 


is  displayed  in  Figure  2  only  for  the  standard  linear  solid.  Note  that  g 
is  singular  at  Y=1  just  as  the  derivation  of  (2)  is  valid  only  for 
0  <_Y  <  1.  A  final  remark  along  these  lines  is  that  if  a  simple  power-law 

behavior  in  time,  i.e.,  g(t)  =  (t/i)  n,  is  assumed,  then 

g(v/a^)  =  ctnr(l-n).  Thus  the  simpler  power-law  model  provides  good 
agreement  with  the  more  realistic  power-law  model  in  predicting  the  ERR, 
yet  has  the  desirable  feature  of  providing  a  model  which  leads  to  a 
boundary  value  problem  in  which  field  quantities  are  more  easily  calculated 
[1]. 

This  paper  concludes  by  addressing  the  effect  of  a  small  perturbation 
in  the  spatial  inhomogeneity  upon  the  ERR.  This  analysis  is  most  easily 
carried  out  in  terms  of  displacements.  If  <u,v>  denotes  the  displacement 
field  corresponding  to  the  original  boundary  value  problem  discussed  above, 
then  by  applying  a  Fourier  transform  in  the  variable  x  to  the 
stress-strain  laws,  the  equations  of  equilibrium  and  the  boundary 
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conditions,  the  boundary  value  problem  corresponding  to  the  modulus 


u(t,y)  =  uc(y/yc) 


g(t) 


may  be  written  as 


L [u,v] (p,y)  +  y  B [u,v](p,y)  =  0 

lim(y/y  )y  B  [u,v](p,y)  =  0 
y^O 


y>0 ,  -00  <  p  <  °° 

-00  <  p  <  00 


where 


lim(y/y  )  B  [u,v](p,y)  =  f(p)  p<0 

y+0 


v(p,0)  =  0 


p>0 


L[u,v]  = 


2  a 
3  u 


.2- 
o  v 


ip  3v  2(l-v)p  u 

1-2V  7^  l-2v 


ip  3u  ( l-2v )p^  v 


2  2(  1-v )  7y  2(l-v) 


Bi^u,v^  =  jy  (p*y)  "  ip  v(p»y) 


B?[u,v]  =  I7  (p,y)  -  ..ipV.  u(p,y) 
y  (1-V) 


f(p)  = 


(l-2v)  o"(p) 
(1-v)  2u c  g(ivp) 


and 


Blu.v] 


Bj (u,v) 
B2(u,v) 


(24) 
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If  the  spatial  inhomogeneity  is  modified  as  e  +  (y/y  ),  the  above 


boundary  value  problem  becomes 


(y/yc)Y  [l[u,v]  +  (Y/y)  8[u,v]]  +  e  L[u,v]  =  0 

•  lim(e  +  (y/y  )YJ  B  [u,v]  =0  -°°  <  p  <  00  (25) 

y+0  1 

lim  (e  +  (y/y  )Y]  B  [u , v ]  =  f(p)  p<0 
y+0  c  z 

w  u( p ,0)  =  0  p>0. 

If  the  solution  of  (25)  is  expanded  in  the  asymptotic  series 


<u,v>  =  l  £n<u  ,  V  >, 
n  n’  n  * 
n=0 


then  it  is  a  straightforward  matter  to  verify  that  substitution  of  this 
series  into  (25)  shows  that  <Ug ,  Vg>  solves  the  unperturbed  problem  (24) 
and  hence  corresponds  to  the  solution  obtained  in  [1].  Having  found 

<u0’  v0'>’  one  solves  for  <un>  vn>  recursively.  For  example,  <Uj,  Vj> 
satisfies 


(y/yc]Y  Ut^,  v^  +  Y/y  B [u  1 ,  ] ]  +  i[GQ,  vq]  =  0 

li4Bl(u0»  'V  +  ^/yc)Y  Bjtuj,  Vj]]  =0  <  p  < 

lim[B2[u0,  vQ]  +  (y/yc)Y  B2  [ux  ,  v^]  =0  p<0 


u1 (p»0)  =  0 


p>0. 


The  ERR  corresponding  to  the  perturbed  problem  (25)  can  now  be 


expressed  by  the  formal  series 
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G-C  Vv.y  dP 


-  /_»  °fVV0,y  dP  +  £  Too  °fVvi,y  dP  +  ^  /loo  dP  + 


=  Gq  +  0(G )  as  e+0 


where  GQ  corresponds  Co  the  ERR  when  Che  spatial  inhomogeneity  is  given 
Y 

by  (y/y •  Thus  when  the  spatial  inhoraogenei ty  is  characterized  by  a 

modulus  of  the  form  £  +  ( y/y £ )Y  and  e  «  1,  corresponding  to  a 
material  that  offers  little  resistance  to  shear  near  the  crack  plane,  then 

the  corresponding  ERR  agrees  with  the  ERR  predicted  by  the  modulus  (y/yY  ) 
to  an  order  of  e . 
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Legends  for  Figures 


The  normalized  effective  depth,  (y, /a  )Y . 

1  e 

Normalized  ERR  for  the  standard  linear  solid  with 

a  =  .l,  n-io,  5-.01  (-),  .  1( - ),  !(-•-). 


Effective  Depth 
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